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Recently, an interesting class of the direct gauge mediation supersymmetry (SUSY) 
breaking models are proposed, in which the minimum of the potential of the SUSY 
breaking field is determined by the inverted hierarchy mechanism. We consider their 
Q-i, cosmological implications. In this class of models, SUSY breaking field has a very flat 

potential, which may have a cosmological importance. Assuming the initial amplitude 
of the SUSY breaking field to be of the order of the Planck scale, it can be a source 
of a large entropy production. A special attention is paid to the cosmological moduli 
problem, and we will see the cosmological mass density of the moduli field can be 
& ' significantly reduced. 
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1 Introduction 



Low energy supersymmetry (SUSY) has been regarded as one of the attractive new physics 
beyond the standard model, since it may provide a natural explanation to the stability of 
the electroweak scale against radiative corrections. However, contrary to our theoretical 
interests, any superpartner of the standard model particle has not discovered yet, and hence 
SUSY has to be broken in nature. Unfortunately, we do not have a clear picture of the SUSY 
breaking in nature, and the understanding of the origin of the SUSY breaking is one of the 
most important issues in the study of the supersymmetric models. 

In recent years, a new framework for the SUSY breaking, i.e., gauge mediated SUSY 
breaking (GMSB) [Ij], has been attracting many interests, and its phenomenological impli- 
cations have been extensively investigated 0. In particular, in GMSB, SUSY breaking in 
the SUSY standard model (SSM) sector is mediated by the gauge interactions which do not 
distinguish flavors. In this scheme, dangerous off-diagonal elements in the sfermion mass 
matrices are suppressed, and serious SUSY flavor changing neutral current (FCNC) problem 
can be evaded. 

In spite of the phenomenological interests, however, cosmology of GMSB is not fully 
satisfactory. In particular, relic abundances of the gravitino |3|, ffl and the moduli field [H 
have been known to be problematic. These problems are often called "gravitino problem" 
and "cosmological moduli problem." In some sense, they are more serious than the gravity 
mediated SUSY breaking case, and they may be crucial weak points of GMSB. (These issues 
will be reviewed in the next section.) However, one should note that there are rooms to 
solve or improve these difficulties; since these difficulties are based on a kind of "minimal" 
assumption, some of them may be solved or relaxed by a new idea. For example, thermal 
inflation |J is proposed to dilute the unwanted particles. 

In this paper, we would like to propose a new mechanism for a large entropy production, 
which can be a resolution to the cosmological problems in GMSB. Our scenario is based 
on a class of models with direct gauge mediation in which the messenger particles have a 
direct coupling to the original SUSY breaking field. In particular, recently, several direct 
gauge mediation models are proposed in which SUSY breaking field has exactly flat potential 
at the tree level 0, |], [| [Kj. (For other classes of models of direct gauge mediation, see 



Refs. [11, [nj.) In this class of models, minimum of the potential of the SUSY breaking 
field is determined by the inverted hierarchy mechanism [jl3| , and SUSY breaking field has 
a very flat potential even after the potential is lifted. In this case, SUSY breaking field may 
play an important role to dilute unwanted particles; with an assumption that the SUSY 
breaking field has an initial amplitude of the order of the Planck scale, various cosmological 
problems can be naturally solved. One virtue of this scenario is that the source of the large 
entropy production is already in the framework of the SUSY breaking mechanism. Therefore, 
the scenario is fairly economical, and we do not have to introduce any new field only for 
the entropy production (like "flaton"), contrary to the case of thermal inflation ||. In the 
following sections, we see how this works, and consider if we may have a cosmologically 
consistent scenario. 
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The organization of this paper is as follows. In Section 0, we give an overview of the 
cosmological difficulties in GMSB. Then, in Section || we briefly review the important aspects 
of the direct gauge mediation model with the inverted hierarchy mechanism. In particular, 
the SUSY breaking field plays a very important role in our discussion, so we see the properties 
of the potential of the SUSY breaking field in some detail. In Section [|, cosmology based 
on the direct gauge mediation model with the inverted hierarchy mechanism is discussed. In 
particular, we concentrate on the cosmological evolution of the SUSY breaking field, and we 
discuss how it improves the cosmological difficulties. Section ^| is devoted to discussion. 

2 Overview of the Cosmology of GMSB 

Before discussing the cosmology of the direct gauge mediation model, let us first briefly 
overview the cosmology of the gauge mediation model. 

Cosmologically, one important outcome of the GMSB is the light gravitino; the gravitino 
mass in this scheme has to be much lighter than the SSM scale mssM- This is because 
mass squared matrices of sfermions have off-diagonal elements of 0(my 2 ). If this effect is 
comparable to the dominant contribution from GMSB, dangerous SUSY FCNC problem 
arises again, which spoils the important motivation of GMSB. Consequently, the gravitino 
becomes the lightest superparticle (LSP) in this scenario. 

Keeping this feature in mind, one of the most famous cosmological constraint in GMSB 
is from the mass density of the gravitino in the Universe. If the gravitino is thermalized 
in the early Universe, and if it is not diluted, its mass density may significantly contribute 
to the energy density of the Universe. Without dilution, mass density of the gravitino is 
proportional to 7713/2, and the Universe is overdosed if the gravitino mass is heavier than 
about 1 keV 0. 

With an enough entropy production after the decoupling of the gravitino from the ther- 
mal bath, gravitino mass heavier than 1 keV may be also viable. However, even in this 
case, the gravitinos are produced in the thermal bath due to scattering and decay processes. 
These secondary gravitinos also contribute to the mass density of the Universe, and hence 
gravitino production after the entropy production has to be inefficient. Since the gravitino 
production is more effective for higher temperature, we obtain an upper bound on the max- 
imal temperature of the Universe after the late entropy production in order not to overdose 
the Universe [f|]. Notice that the interaction of the longitudinal gravitino with matter is 
proportional to Tn^, 2 , and hence the constraint becomes weaker for a heavier gravitino. In 
Fig. |l|, we show the upper bound on the maximal temperature T max as a function of the grav- 
itino mass; the upper bound is from ~ 100 GeV to ~ 10 8 GeV, for the gravitino mass 1 keV 
- 1 GeV. In particular, the constraint is very strict in the case of 1 keV ^ m.3/2 ^5 100 keV, 
where the gravitinos are mainly produced by the decay processes. If the Universe starts with 
a temperature higher than this upper bound, the Universe is overdosed by the gravitino, 
unless there is an enough entropy production below this temperature. 

More serious problem arises in the framework of superstring models. In superstring 
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Figure V. Upper bound on the maximal temperature of the Universe as a function of the 
gravitino mass 7723/2- 

models, dilaton and moduli fields (which we call "moduli" fields hereafter) exist which are 
the flat directions related to symmetries in the superstring theory. Mass of the moduli 
field originates to SUSY breaking effect, and is of the same order of the gravitino mass. 
(Throughout this paper, we approximate the mass of the moduli field to be the gravitino 
mass 7773/2.) Generically, moduli field takes an initial amplitude of the order of the Planck 
scale, unless our vacuum lies at or near a point of enhanced symmetry [0.0 Then, it starts 
oscillation when the expansion rate of the Universe becomes comparable to the mass of the 
moduli field. Since the interactions of the moduli field are suppressed by inverse powers of 
the Planck scale, moduli field lighter than about 100 MeV has a lifetime longer than the 
present age of the Universe. In GMSB, this is (almost) always the case. In this case, mass 
density of the moduli field becomes enormous, if there is no dilution. Assuming the radiation 
dominance before the moduli field starts to move, naive calculation results in the density 
parameter of the moduli field as 

~ 6 x 10" x f^r /4 p^v /2 ( m 2 , (2 .D 

UOOy VIOOkeV/ \M m J ' V ; 

where h is the Hubble constant in units of 100 km/sec/Mpc, g* is the effective number of 
the massless degrees of freedom when the moduli field starts to move, M* ~ 2.4 x 10 18 GeV 
is the reduced Planck scale, and <J)q is the initial amplitude of the moduli field. For example, 
for 7773/2 = 100 keV, the initial amplitude has to be smaller than ~ 10 _7 M* in order not 

^However, it is difficult to construct a realistic model which has our vacuum as a symmetry enhanced 
point. 
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to overdose the Universe, and this is an extreme fine tuning of the initial condition. In 
other words, if the initial amplitude takes the natural value (i.e., ~ M*), a large entropy 
production is inevitable for a viable cosmological scenario. 

Even if we adopt a large entropy production to dilute the gravitino and the moduli field 
away, it is still non-trivial whether a consistent cosmological scenario can be obtained. If 
there is a large entropy production, it also dilutes the possible baryon asymmetry generated 
in the early stage of the Universe. Naively speaking, baryon number asymmetry has to be 
generated after the entropy production. However, in some case, the reheating temperature 
after the entropy production becomes too low for baryogenesis. In Ref. ||, it has been 
pointed out that Affleck-Dine mechanism for baryogenesis [JHJ may be able to generate 
enough baryon asymmetry even if there is a large entropy production. This topics is reviewed 
later. Candidates of the cold dark matter (CDM) is another interesting issue in GMSB. In 



the gravity mediated SUSY breaking scenario, the lightest neutralino [I7j or sneutrino |T 
can be a promising candidate of the CDM, if it is the LSP. However, in GMSB, they cannot 
be the CDM, since they can decay into gravitino and their superpartner. Several candidates 



of the CDM are proposed in the framework of GMSB |19| , |20 |, but these candidates are also 
diluted by the entropy production. One interesting candidate is the coherent oscillation of 
the moduli field, if its energy density is diluted enough. Therefore, it is very important to 
consider the possibility to dilute the energy density of the moduli field down to fi^ ^ 1. 

Another class of cosmological problems are related to the structure of the scalar potential; 
the scalar potential may have unwanted minimum which is deeper than the phenomenolog- 
ically viable local minimum. For example, original low energy gauge mediation model may 
have a color breaking minimum pl| , and the models proposed in Refs. [§, ^ [K| have a SUSY 
preserving true vacuum at the origin of the potential of the SUSY breaking fieldQ Usually, 
the tunneling rate to the true vacuum can be so small that the transition does not happen 
for the time scale of the age of the Universe. Therefore, it is phenomenologically consistent 
once the SUSY breaking field is trapped in the minimum we want. However, cosmologically, 
we have not understood how the SUSY breaking field is trapped in the relevant (local) min- 
imum, not in the unwanted (global) one. In particular, if we assume a naive SUSY breaking 
phase transition, many horizons choose the unwanted deeper minimum, and the current 
horizon contains many regions which dropped into the unwanted minimum. Notice that it 
is unclear whether the thermal inflation could solve this problem, even though the reheating 
temperature after the thermal inflation is relatively low. This is because the current horizon 
scale contains many different horizons before the thermal inflation. Therefore, even if the 
SUSY breaking phase transition occurs before the thermal inflation, current horizon still 
contains many regions of unwanted minimum. 

Keeping these arguments in mind, it is important to develop a cosmologically consistent 
scenario based on GMSB. Importantly, one should remember that above problems are usually 
based on "minimal" assumption, and in particular in direct gauge mediation model, the above 
arguments do not take account of a possible effect from the SUSY breaking field. In the 



# 2 However, the structure of the potential is model-dependent, and models without unwanted minimum 
may be constructed. 
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following sections, we will see what happens if we include its effects. 

3 Model 

In this section, we first briefly review the class of models we are interested in. As we 
mentioned in the introduction, we consider a cosmology of the direct gauge mediation model 
in which the potential of the SUSY breaking field is stabilized by the inverted hierarchy 
mechanism. In this section, we discuss the general features of such models |7], ||, §, |K|. An 
explicit example of the model is shown in Appendix [A|. 

3 . 1 Framework 

The model is based on the symmetry G = Gs x Gb x Gsm- Here, Gs is the strong gauge 
interaction whose dynamics induces the gaugino condensation. On the other hand, Gb is 
introduced to stabilize the minimum of the potential of the SUSY breaking field, and Gsm 
is the standard SU(3)c x SU(2)l x U(1)y gauge croup. In this framework, we introduce the 
following class of chiral superfields in order to induce a desired dynamics to break SUSY: 
the SUSY breaking field £ which transforms under Gb, "messengers" q + q which transform 
under Gb and Gsm, and "strongly interacting" fields Q + Q which transform under Gs and 
Gb- These chiral superfields have a superpotential of the form 

W = y Q ZQQ + y q Zqq. (3.1) 

Notice that this superpotential has i?-symmetry which is non-anomalous for Gs; under this 
-R-symmetry, E has charge +2. As we will discuss later, this i?-symmetry is explicitly broken 
by supergravity effect. 

We construct a model so that S has a flat direction which is parameterized by an invariant 
made out of E. On this flat direction, Gb is broken by the vacuum expectation value (VEV) 
of S. The first requirement to the model is that, on this flat direction, all but one degrees 
of freedom in £ are eaten by Higgs mechanism. We parametrize the remaining degrees of 
freedom by X. 

On this flat direction, chiral superfields coupled to S also acquire masses. This effect 
changes the /3-function of Gs at fi ~ yq(E), since Q and Q have mass of tuq ~ 2/q(E) and 
decouple at this scale. As a result, the strong scale for the theory above tuq, A, differs 

from that after Q and Q decouples, A e g. These two quantities are related as A^ Gs = 
A 3/iG s ra^Q , where //q s and \iq are Dynkin indices for the adjoint and Q+Q representations 
of Gs, respectively. The second requirement to the model is that these Dynkin indices 
satisfy the relation fiQ s = [Mq, so that the linear superpotential is induced by the gaugino 
condensation; below the strong scale, effective superpotential is induced by the gaugino 
condensation as 

Wes ~ Ae ff ~ {ygXyi 1 ^ a3-mq/mg s _ (3.2) 
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Therefore, if /xg s = hq, superpotential is linear in X, and supersymmetry is broken by the 
VEV of the F-term of the X field; Fx ~ 2/qA 2 . (This class of model of SUSY breaking is 
originally discussed in Ref. fll4| .) 

With the above superpotential, scalar potential is given by 

v _ \dxW eS \\ : y 2 Q ^ 

d x d x K ~ zv(x*,xy { • ] 

where dx represents the derivative with respect to X, K is the Kahler potential, and Z^ 
is the wave function renormalization of S. At the tree level, Z^ = 1, and hence V does 
not depend on X. In this case, VEV of X is undetermined. However, once we include the 
radiative corrections, the situation changes. Since £ interacts through gauge and Yukawa 
interactions, non-trivial Kahler potential is induced by radiative correction. X dependence 
of Z-% is governed by the renormalization group equation (RGE). At the one loop level, RGE 
for Zs is given by 



dlnZ-s 1 



(iln/x 167T 2 



C B g 2 M - C Q y 2 Q (ft) - C q y 2 q (fi) , (3.4) 



where coefficients Cb, Cq, and C q are all positive. If the /3-function for Z^ vanishes, we 
have a minimum (or maximum) of the potential. The important point is that the gauge 
and the Yukawa contributions have opposite signs in the /3-function of Zg. As we can see 
in Eqs. (3.3) and (|3.4j ), gauge contribution makes Z^ larger at higher energy and drives \X\ 



to a larger value, while Yukawa contribution affects in the opposite way. As a result, if the 
gauge piece dominates the /3-function in the low energy, and also if the effect of the Yukawa 
terms wins in the high energy, X has a minimum at \X\ = v where the /3-function vanishes: 

C B g 2 B (v)=C Q y 2 Q (v) + C q y 2 (v). (3.5) 

Therefore, in this class of models, potential for X is stabilized by the scale dependence of 
the wave function renormalization factor Z%. At the minimum, F-component of the X field 
has a non- vanishing VEV of 0(?/qA 2 ), and SUSY is broken. Then, fermionic component of 
X becomes the goldstino, and in supergravity, it is absorbed in the gravitino. In this case, 
gravitino mass is related to Fx as 

-s /2 = (3.6) 

Once X gets a VEV v, q and q acquire a SUSY preserving mass of ~ y q v, as well as 
SUSY breaking mass squared for the scalar component of ~ y q Fx- Since q and q have 
quantum numbers under the standard model gauge group, the SUSY breaking masses for 
the SSM superparticles arise by integrating out these messenger particles. As in the case of 
well-known ordinary gauge mediation model, the ratio, 

B Q = ^, (3.7) 
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Figure 2: Contours of the constant gravitino mass on the v vs. Bq plane. 



determines the scale of the SUSY breaking masses in the SSM superparticles which are of 
0(«sm/ '4vt)_Bq, with a SM being the appropriate gauge coupling of the standard model gauge 
group. Bq should be in the range of 10 4 GeV < B Q < 10 5 GeV; the lower bound is from 
experimental lower bounds on the masses of the superparticles, while the upper bound is 
from the naturalness point of view. 

In the following discussions, we use v and Bq as independent parameters in the model, 
and rewrite other quantities as functions of them. For example, the gravitino mass is given 
by 

In Fig. [|, we show the contours of the constant gravitino mass on the v vs. Bq plane. 

Before closing this subsection, we discuss the allowed range of v. Recalling that ~ y 2 v 2 
and ~ y q Fx are the diagonal and off-diagonal elements of the mass squared matrix of the 
messenger scalars, v has to be larger than ~ Vq^Fj/ 2 for the positivity of the eigenvalues 
of the mass squared matrix; otherwise, messenger scalar has a VEV and the standard model 
gauge groups are broken. Another constraint is from the stability of the SUSY breaking 
minimum. In some class of models, there is a SUSY preserving true vacuum at the origin 
(\X\ = 0). In this case, the tunneling rate from the SUSY breaking vacuum to the true 
vacuum has to be small enough, so that the lifetime of the SUSY breaking vacuum is longer 
than the age of the Universe. This issue is discussed in Ref. ||, and it requires w/A^ 10. 
Since Fx ~ 2/qA 2 , IO^q 1 ^ 2 'f][ 2 is required in models with true vacuum at the origin. 
Notice that this constraint is more stringent than the one from the stability of the messenger 
potential, if y Q ~ y q . 
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One important lower bound on v is derived for the validity of our perturbative approach. 
If v is close to the strong scale A, SU(2)g dynamics generates non-perturbative Kahler 
potential, and our perturbative arguments break down. At the scale A e fj, we expect a non- 
perturbative contribution to the Kahler potential of the form 5K ~ [(W a VV Q ,)(VVaVV a )] 1 / 3 ~ 
A* ff A cff [^, ||, where W a is the spinor superfield for SU(2) S . With the naive dimensional 
analysis p3| , coefficient of this operator is estimated to be of 0((47r)~ 2 / 3 ). Requiring this 
effect to be smaller than the perturbative one, v has to be larger than 10 9 — 10 10 GeV ||. 
However, it may be possible that the minimum of the potential exists even with the non- 
perturbative effect. (For a model in which the SUSY breaking minimum is stabilized by 
non-perturbative effects, see Ref. In this case, we need to make some dynamical 

assumptions not to change our following discussions, since the non-perturbative effect is not 
calculable. However, most importantly, strong dynamics does not break the i?-symmetry, 
and hence the properties of the i?-axion is unchanged. Consequently, evolution of the R- 
axion is basically unaffected, even if the non-perturbative effect becomes important at the 
minimum of the potential. With these caveats in mind, we also consider the case with 
v ^5 10 9 GeV in the following discussion. In this case, we make an implicit assumption that 
the minimum of the potential is stabilized even after taking account of the non-perturbative 
effect. 

Finally, we comment on the upper bound on v. As can be seen from Eq. Q3.8p , gravitino 
mass becomes larger than ~ 100 GeV for v ^ 10 16 GeV with Bq ~ 10 4 GeV. In this case, 
supergravity contribution to the sfermion masses becomes comparable to the gauge mediation 
piece, and hence SUSY FCNC problem arises again. Therefore, v <C 10 16 GeV is required 
to solve SUSY FCNC problem. 



3.2 Potential of X 

In order to discuss the cosmological implication of X, it is important to understand the 
properties of its potential, which is the subject of this subsection. 

We are interested in the models in which the minimum of the X field is determined by 
the inverted hierarchy mechanism. In this class of models, the potential at the global level 
is only lifted by the renormalization group effects, and hence it can be expanded by powers 
of \nX*Xf^\ 



Vo 



global 



F 2 X /Z 1 



F 2 



167T 2 

d(M) ln ^ + 



16tt 2 




z 2 M 




x*x 



+ 



(3.9) 



# 3 After the second line, we omit the constant piece. We assume that the cosmological constant is cancelled 
out by supergravity effect, and hence the constant term does not matter. 
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with 



2 dZ% _ 2 2 d Zy, 



Z'x = 8vr 2 -^, Zl = (Stt 2 ) 2 -^^. (3.10) 

Here, ( n is a 2n-th polynomial of gauge and Yukawa coupling constants. We factorized rele- 
vant powers of lQir 2 , so that coefficients in ( n become typically of 0(1). Without cancellation, 
( n is close to 1 if some of the gauge or Yukawa coupling is close to 1. 

In Eq. (|3.9| ), fi can be an arbitrary scale; \i dependence is cancelled out by the renormal- 
ization group effect. If we take fi = v where Z^ vanishes, V^iobai starts with a term which is 
quadratic in lnX*X, and hence the potential has an extremum at |Jf| = v. In particular, if 
(2(v) > 0, potential has a minimum there, which is what we desired. Hereafter, we assume 
that the gauge and Yukawa coupling constants are arranged so that (2{v) > 0. Around this 
minimum, potential starts with (lnX*X) 2 term which is suppressed by (167r 2 )~ 2 . However, 
once |X| becomes much larger (or smaller) than v, Ci at that scale may become large so 
that the potential is approximated to be linear in lnX*X, with being suppressed only by 

(167T 2 )- 1 . 

So far, we have discussed the potential in the framework of global SUSY. However, 
supergravity effect also generates important terms in the potential. First of all, all the 
scalar fields receive SUSY breaking masses of the order of the gravitino mass. This effect 
becomes important especially when the amplitude of X becomes large. Another important 
implication is that the i?-symmetry is (explicitly) broken due to the supergravity effect if 
the cosmological constant is cancelled out pifl . Under U(1)r symmetry, superpotential has 
the charge of +2. However, in order to cancel the cosmological constant, a constant term 
is needed in the superpotential. From the cross term between them, /^-symmetry breaking 
potential is induced: 

Vr ~ "||(** + X) x /(X*A7M 2 ), (3.11) 
where / is an unknown function. We expand this function as 

f(x) = k + hx + ---, (3.12) 

where coefficients k n are expected to be of 0(1). By combining these contributions, super- 
gravity contribution is of the form: 

^sugra ~ m 2 3/2 X*X + Vr. (3.13) 

In fact, the linear term in Vr may cause a problem. If we neglect the logarithmic terms, 
non- vanishing k shifts the minimum of the potential from X = to X ~ k^M*. Thanks to 
the logarithmic term in V^iobai, potential can have a minimum at \X\ = v <C M*. However, 
with the potential 

V ~ ^F 2 X \nX*X - k j^(X + X*) + mj /2 X*X + • • ■ , (3.14) 
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another minimum still exists at |X| ~ M* when ko ~ 0(1). With this minimum, X is more 
likely to settle down to this unwanted minimum if X has an initial amplitude of 0(M*). 
This is because the potential is dominated by the supergravity contribution for such a large 
amplitude, and hence X does not feel the effect of the logarithmic piece unless |X| becomes 
small enough. In order to remove this unwanted minimum, we assume k to be small enough, 

fc <|. (3.15) 

If d is of 0(1), this is a tuning of 10 % level, and we believe accidental cancellation would 
be enough for this suppression.^ 



3.3 Properties of the Physical Modes 

At around the minimum (\X\ ~ v), we have two physical scalars from X. In order to discuss 
the properties of these fields, it is convenient to parametrize the X field as 



X=[v + -=a )e ia/Vlv . (3.16) 




Expanding the potential around the minimum, we obtain the mass of the a as 

2 _ Z'Hv) F 2 ( 2 (v) Fl 



(16tt 2 ) 2 v 2 (16tt 2 ) 2 v 2 

where we normalized as Z%(v) = 1. Notice that a is as heavy as the SSM superpartners if all 
the gauge and Yukawa coupling constants are of the same order. On the contrary, a is the 
pseudo-Nambu-Goldstone boson for the i?-symmetry, which is usually called i?-axion. The 
main source of the i?-axion mass is the i?-symmetry breaking term from the supergravity 
effect; 

V ft = -k jt(X* + X), (3.18) 



where k is (unknown) 0(1) constant introduced in Eqs. (|3.11|) and (|3.12|) . (Around the 



minimum \X\ = v, higher order terms are suppressed by powers of v 2 /M 2 .) From this 
potential, the i?-axion mass is given by [p4|p°| 



.2 _ ;. F X 



m a = k 



6 GeV x k, 



1/2 



Bn \ ( v \V2 



10 5 GeV/ V10 10 GeV 



(3.19) 



# 4 The values of ko at \X\ ~ M* and at |X| ~ v are supposed to be different. In particular, in our case, 
moduli field exists which has a very large initial amplitude of 0(M*). The value of ko should be affected by 
the evolution of the moduli field. 

^ 5 There is another contribution to the i?-axion mass from the QCD anomaly, which is, however, much 
smaller than the supergravity effect. This fact suggests that it is difficult to use this i?-axion as a solution 
to the strong CP problem, unless the supergravity effect is much smaller than the naive expectation. 
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(3.20) 



Next, we consider the decay rate of these fields. The i?-axion dominantly decays into 
gauge boson pairs.|^] Since the i?-axion couples to the messenger multiplets q and q which 
transform under SU(3)c x SU(2) L x U(l)y, the i?-axion has a coupling to the standard 
model gauge bosons at the one-loop level. Then, its decay rate is calculated as 

r _ n G (asubXml 
a 32tt V 4tt J v 2 ' 

Here, uq is the number of the final state, b q is the /3-function coefficient of q + q, and «sm 
is the corresponding gauge coupling constant. For example, for the decay into the gluon 
pair (a — > gg), ric = 8, b q = N$, and «sm = ot s , while for the decay into the photon pair 
(a — > 77), n G — 1, b q — |iV5, and «sm = « em , with N 5 being the number of messenger chiral 
superfields in units of 5 + 5 representations of SU(5). (For the model shown in Appendix |A|, 
N 5 = 2.) 

On the other hand, the a field dominantly decays into the i?-axion pair (a — > aa) . From 
the following Lagrangian: 

£ = d^X*d^X ~ -d^ad^a + ^d^ad^a + -^-ad^ad^a H , (3.21) 

we obtain the decay rate of this process to be 

1 

64tt v 2 

Notice that a also decays into the gauge boson pairs, and the decay rate for this process 
is given by a similar formula as Eq. (|3.20| ) with m a being replaced by m a . Comparing the 



r. = — -f ■ (3-22) 



decay rates for these processes, we can see that a —> aa is the dominant decay mode for a. 

The decay mode into the gravitino pair is potentially significant. Indeed, the SUSY 
breaking field has an interaction as 

with ipp being gravitino. This interaction induces the decay process X — > ip^ip^, and in 
particular, the decay into the longitudinal component is the most important. With the 
Lagrangian given above, we obtain the decay rate as 



rn 



5T X = T(X - ifiM = — , (3.24) 



32tt F x 



# 6 If the ordinary quarks and leptons have non-vanishing i?-charge, the i?-axion may decay into these 
fermions. However, it is unknown whether the i?-symmetry can be consistently defined in the SSM sector, 
and the couplings of the i?-axion to the quarks and leptons are model-dependent. In particular, these 
couplings cannot be fixed unless we specify the mechanism to generate so-called fi- and B-parameters. 
Furthermore, these decay processes are chirality suppressed. In this paper, we assume that these decay 
modes are negligibly small. The i?-axion may also decay into gaugino pairs. However, the i?-axion is lighter 
than the gauginos in most of the parameter region we consider. Therefore, we do not consider this decay 
mode. Notice that the decay rate for this process is at most comparable to that into gauge bosons. Therefore, 
qualitatively, the following arguments are unchanged even if the i?-axion decays into gaugino pairs. 
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with X = a and a. Comparing Eq. ( |3.24| ) with Eqs. ( |3.20| ) and Q3.22| ), the branching ratio 
for the process X — > ip^ip^ is estimated as 

Br(a -> ~ ^ ~ J- f ? ™' ) 2 , (3.25) 

T a n G \(a SM b q /4:7r)BQ J 

Br(a^^)^^2^\ (3.26) 

Notice that these branching ratios are much smaller than 1. In the parameter region we 
consider, the i?-axion mass is usually smaller than 0(«sm/ '4k)Bq. Furthermore, as shown 
in Eq. ( |3.17 ), m CT is smaller than Bq. Therefore, in both cases, the decay modes into the 



gravitino pair are suppressed. Even with small branching ratios, however, the decay of the 
.R-axion and a may overproduce the gravitino, resulting in too much mass density of the 
Universe. We will come back to this point in the next section. 



4 Evolution of the Scalar Fields 

Now, we are ready to discuss the evolutions of the scalar fields, i.e., the SUSY breaking field 
X and the moduli field 0. Since the behavior of X changes at \X\ ~ v, we first consider the 
evolution when \X I ^> v. Then, the behavior of X when \X\ ~ v is considered. 



4.1 \X\ ~ M. 



* 



In our discussion, we adopt the picture of the inflationary cosmology; we assume a primordial 
inflation which solves horizon and flatness problems. During this inflation, X and are 
shifted from their minimum, and these fields have very large initial amplitudes. We assume 
that their initial amplitudes are both of 0(M*). This may happen because of the chaotic 



assumption on the initial condition of scalar fields | 25| , due to modification of the scalar 



potential during the inflation with large expansion rate or in the case of no-scale type 



supergravity [2f| . Importantly, e-folding of this primordial inflation is large enough to solve 



the horizon problem, and hence X and have coherent initial states for the scale of the 
current horizon. After the primordial inflation, reheating happens, and the Universe becomes 
radiation dominated. At this stage, X and keep their initial values as far as the expansion 
rate is larger than their masses. 

Once the expansion rate of the Universe becomes comparable to the mass of X and 
0, these fields start to move. Since these fields (in particular, X) are assumed to have 
large amplitudes of 0(M*), the potential for these fields are initially dominated by the 
supergravity contributions. Furthermore, with these large initial amplitudes, energy density 
of the radiation becomes comparable to those of X and when these fields start to move. 
Then, after this stage, energy density of radiation decreases faster than those of X and 0, 
and the Universe is dominated by the coherent oscillation. 
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Once the moduli field starts oscillation, its evolution is quite simple. Approximating 

3/2 ( 



its potential to be quadratic as ~ ml /2 (f) 2 , the energy density of the moduli field scales as 



R~ 3 , with R being the scale factor. 

One may worry about a possible shift of the minimum of the moduli potential, in par- 
ticular, due to the non-vanishing expansion rate H induced by X. When the Universe is 
dominated by a condensation of a scalar field (in our case, X), supergravity effect induces 
extra terms in the moduli potential, which are proportional to H 2 [15]. This effect shifts 



the minimum of the moduli potential, and <p oscillates around the shifted minimum. If the 
minimum moves more drastically than 0, behavior of (f) is affected by the change of the 
potential, and our argument breaks down. However, in our case, this is not the case; time 
scale of the shift of the minimum is H~ l , which is much longer than that of the oscillation 
JTiT^. As a result, can catch up with the shift of the minimum. Therefore, if we consider 
the oscillation around the shifted minimum, our argument is unchanged. In particular, this 
does not affect the calculation of the current energy density of the moduli field since the 
shifted minimum smoothly approaches to the true minimum as the Universe expands. (For 
details, see Appendix [B].) 

Evolution of the SUSY breaking field X is more complicated. When X starts to move, 
angular component of X (i?-axion mode) is excited by the Affleck-Dine mechanism and R- 
number is generated. Furthermore, when logarithmic potential takes over the supergravity 
contribution, energy density of X decreases much slower than that of the moduli field. 

Let us first discuss the generation of the .R-number due to the Affleck-Dine mecha- 
nism | 16fl . For our argument, it is convenient to define the .R-number density: 

n R = i(X*X - X*X). (4.1) 

Then, with the .R-symmetry breaking terms given in Eq. (|3.11| ), time evolution of n R is given 
by 

( dV dV \ 

nR+3Hna = _^_ A .__ A ..j 

= 2F 2 x lm(X/M.) x f(\X\ 2 /M 2 ,). (4.2) 



Notice that, if the potential respects the .R-symmetry, the right-hand side of Eq. ( |4.2| ) van- 
ishes. 

In our case, .R-symmetry breaking effect is the largest when X has the maximum am- 
plitude, and .R-number asymmetry is generated when X starts to move (see Appendix [C]) . 
.R-number at this stage is estimated as 

n R ~ H- 1 x 2F 2 lm(X /M*) x f(\X \ 2 /M 2 ). (4.3) 

By using H ~ m 3 / 2 , we obtain 

2^sinfl 

n R , (4.4) 

m 3/2 
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Figure 3: Initial motion of the SUSY breaking field X on the complex X plane. We choose 
the set of parameters as X = M*e 47r//2 , £i = 0.3, k\ = 0.2, and other C n 's and fc n 's are taken 
to be 0. 

where the initial value of X is parameterized as Xq = \Xo\e ie °, and £ ~ (|Xo|/M*) x 
/(|X | 2 /M2) is expected to be of 0(1) if \X \ ~ M*. 

In Fig. |3], we show a typical behavior of X on the complex X plane. As an example, we 
choose the set of model parameters as 

X = M,e in/ \ Ci = 0.3, fc! = 0.2, (4.5) 

and other £ n 's and A; n 's are taken to be 0. 

In order to discuss the evolutions of the moduli field and X simultaneously, it is convenient 
to take the ratio of the energy density of <fi, p$, to Ur. Since the initial energy density of the 
moduli field is of the order of ml^o ~ Fx x (^o/-^*) 2 , we obtain 

n R 2£sin^ \Mj ' 1 ' J 

Notice that and are both proportional to -R -3 , and the ratio p^/nn remains constant 
as far as the i?-symmetry breaking effects can be neglected. 

4.2 v<\X\<M, 

Once the SUSY breaking field X and the moduli field start to oscillate, their amplitudes 
adiabatically decrease with the expansion of the Universe. Their relation is determined by 
Eq. (|4.6|), and the final result can be derived without discussing the detail of their evolutions. 
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However, in this subsection, we discuss how their amplitudes behave when v^z \X\ ^ M* for 
a better understanding of the behaviors of X and 0. 

As discussed in Appendix [D], amplitude of the scalar field (p (corresponding to X and 0) 
obeys the relation 

m e s(p 2 R 3 = const., (4.7) 
where the "effective mass" m e g is defined as 

2 1 / 1 dV 1 dV \ . , 

Effective masses of X and depend on their amplitudes differently. Potential for the 
moduli field is parabolic, and hence m e s is independent of the amplitude. Thus, the 
amplitude of the moduli field scales as R~ 3 / 2 . 

We next discuss the evolution of X. For this purpose, let us remind the structure of the 
potential of Xip 7 ] 

V ^Fl In X*X + m\ l0 X*X. (4.9) 

167T 2 * 3/2 V ) 

As mentioned in the previous subsection, supergravity contribution wins the global SUSY 
contribution when \X\ is large. Comparing two contributions, supergravity effect is more 
important above a threshold value X t h r which is estimated as 



X thl ~ ^M*. (4.10) 

For \X\ ^ X t hi, supergravity effect wins the global SUSY contribution, and m e fj ~ m 3 / 2 for 
X. In this case, the amplitude of X scales as R~ 3 ^ 2 . On the other hand, for \X\ ^ X t hr, 
potential is dominated by the logarithmic piece, and m e fj oc (Al" 1 . In this case, the amplitude 
of X is proportional to R~ 3 ^_ 



Comparing the evolutions of X and 0, their amplitudes are related as 

a/\y\ / ^o/|A | : |X|>X thr (Ait's 

1 (A thr /I^l) 1/2 x (0o/|X o |) : \X\<X thl ' 

with 0o an d X being the initial values of X and 0, respectively. For \X\ ^X t h r , X and 
scales in the same way, while for \X \ ^ Xthi, the amplitude of X decreases faster than that 
of 0. 

The evolution of the energy density of X and is also non-trivial. When \X\ ^ X t h r , both 
X and obey parabolic potential, and their energy density scale as R~ 3 . However, once 



# 7 We neglect the i?-symmetry breaking part; even if we include its effect, the result is unchanged. 
# 8 If \X\ becomes close to v, potential is approximately proportional to (lnX*X) 2 . In this case, this relation 
receives a logarithmic correction; \X | (In IXI) 1 / 2 ^? 3 = const. 
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the amplitude of X becomes smaller than X t hr, potential for X is lifted only logarithmically, 
and energy density of X decreases very slowly with the decrease of \X\. The important 
point is that, once the potential of X is dominated by the logarithmic piece, energy density 
of (f) decreases much faster than that of X, and hence the energy density of the Universe is 
dominated by X. Therefore, when X decays, there can be a large entropy production to 
dilute the moduli field away. 

When the amplitude of X becomes of 0(v), X is trapped in the minimum of the potential. 
Then, the amplitude of X is fixed to be v, and we need to consider the excitations around the 
minimum. They are phase degrees of freedom (i.e., i?-axion) and radial degrees of freedom. 
Effects of these fields are discussed in the following subsections. 

Some of the models have a SUSY preserving true vacuum at the origin X = 0, and one 
may worry whether X can be smoothly trapped in the SUSY breaking (false) vacuum of 
\X\ = v. If X would overshoot down to the origin, this scenario would not be phenomeno- 
logically viable. Since X follows an elliptic trajectory on the complex X plane, as shown in 
Fig. |3|, X can result in the SUSY breaking minimum at \X\ — v without being affected by 
the potential for \X\ <C v. Important point is that, thanks to the .R-number conservation 
in the comoving volume, the orbit of X is always elliptic once the i?-number is generated. 
Consequently, even though the amplitude of X decreases with the expansion of the Universe, 
orbit does not pass by the origin if enough .R-number asymmetry is generated. In this case, 
when \X\ ~ v, X traces a trajectory which is (approximately) parallel to the minimum of the 
potential (\X\ = v), and X is smoothly trapped in the SUSY breaking vacuum. Therefore, 
once enough Renumber asymmetry is generated, X is expected to result in the minimum of 
the potential |X| = v, irrespective of the structure of the potential for \X\ <C v. Of course, 
this scenario depends on the initial value of the Renumber asymmetry. Importantly, if the 
initial amplitude of X is as large as M*, Affleck-Dine mechanism can generate very large 
asymmetry, as discussed in the previous subsection. With a reasonable choice of the model 
parameters, we can easily have an elliptic trajectory of X with eccentricity close to 0. This 
fact helps us to understand how X can be trapped in the SUSY breaking vacuum at \X\ — v. 

4.3 Effect of the i?-axion 

Once the amplitude of X becomes of 0(v), the X field is trapped in the minimum \X\ = v. 
After this stage, the imaginary part (i?-axion) and the real part of X behave differently. 
Thus, we discuss their effects separately. 

First, we argue the effect of the i?-axion a. Since a is the pseudo-Nambu-Goldstone 
boson, its potential is approximately flat with slight perturbation due to V#. When we can 
neglect the i?-symmetry breaking effects, motion of a corresponds to the phase rotation of 
X. At this period, it is convenient to parametrize the X field as 

X = ve iuJt . (4.12) 

Here, we neglect the real part of the X field. Then, the Renumber is given by 

n R = 2uv 2 = 2p a /u, (4.13) 
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where p a is the energy density of the .R-axion. With the expansion of the Universe, ui 
decreases adiabatically; since the .R-number in the comoving volume is conserved, u scales 

as ir 3 .0 

Once the energy density of the .R-axion becomes less than the difference of the potential 
energy AV on the circle \X\ — v, the .R-axion starts to oscillate around its minimum. Once 
this happens, .R-symmetry is not a good symmetry any more, and the .R-axion approximately 
obeys the parabolic potential. After this stage, energy density of a scales as R~ 3 . With the 
potential given in Eq. ( |3.18| ), the difference of the potential energy is given by AV = V(X = 
—v) — V(X = v) = 4m^t> 2 , and hence two cases should be matched when u ~ 0(m a ). 
Connecting two cases at u = u c , we obtain 



P^> P<t> rns/2 



Pa (u c n R /2) (jj c £sm6 




(4.14) 



This ratio remains constant until the .R-axion decays.^] 

Motion of the .R-axion at p a ~ m 2 a v 2 is complicated, and analytic estimation of u c is 
difficult. In our analysis, we numerically followed the evolution of the .R-axion, and estimated 
the value of u for matching. As a result of the numerical calculation, we found that two 
cases should be connected with 

uj = cu c ~ 4m a . (4-15) 

In the following discussion, we use u c = Am a . 

When the expansion rate of the Universe becomes comparable to the decay rate of the 
.R-axion, the .R-axion decays and the Universe is reheated. By using the instantaneous decay 
approximation, the reheating temperature T R is estimated as 



where we used the decay rate given in Eq. ( |3.2U| ) in the second line. (A special case where 



this may be irrelevant will be discussed later.) At this stage, the energy density of the R- 
axion is converted to that of the radiation, and large amount of entropy is produced. At the 

# 9 For a scalar field ip with flat potential, equation of motion is given by (p + ZHip = 0. By solving 
this equation, we obtain ipR 3 = const. Evolution of u is consistent with this relation, if we re-interpret 

^iut y gia/^/2v 

#10 If the amplitude of X is largely fluctuated, domain wall is produced when the i?-axion gets trapped 
in the minimum of the potential. Such a fluctuation is generated during the primordial inflation, and the 
domain wall production may be effective if the expansion rate during the inflation is larger than of 0{v) |2^| . 
Therefore, if we adopt a primordial inflation with relatively small expansion rate, the domain wall production 
can be evaded. Furthermore, even if the domain wall production is effective, collapse of the domain wall 
results in semi-relativistic i?-axions with averaged energy of ~ 3m a p8||. Therefore, domain wall production 



does not modify the ratio p<p/p a given in Eq. (4.14) so much, and the dilution factor calculated below is 
almost unchanged. 



17 



decay time, energy density of radiation is given by p ra( j = p a , and hence 




3 ^R x : ? /2 , 1-^1, (4-17) 



s 4 4m a £sin6> \M* 



where s is the entropy density which is related to the energy density of the radiation as 
Prad — \Ts. We compare Eq. Qi.17 ) with the ratio of the critical density to the entropy 
density in the present Universe: 



Pc 
S 



~ 3.6 x 10~ 9 GeV x K\ (4.18) 



and we obtain 



- 40 x kT (a SM b g n^) f^r 1/4 {^-T ( J f^V , (4-19) 

<$> o V bM « « J VlO/ UO 5 GeV/ UO 10 GeV/ V M */ 

where we used £ sin #0 ~ 1. 

Now, we can calculate the numerical value of the density parameter, and see if it is 
cosmologically viable. For this purpose, we first have to specify the dominant decay mode 
of the i?-axion, since the result depends on the decay width T a . If the i?-axion mass is large 
enough (probably, for m a ^ 1 GeV), the i?-axion decays into the gluon pair. In this case, we 
can use the calculation based on the perturbative QCD, and the decay rate is given in the 
formula in Eq. ( |3.20| ) with «sm = ck s . However, if the i?-axion mass becomes as small as (or 
smaller than) ~ 1 GeV, Eq. (|3.20| ) may not be reliable, since in this case, decay rate into 



multi-meson final states has to be calculated. However, if the i?-axion mass is light enough, 
decay modes into multi-meson final states are kinematically forbidden. Since the i?-axion is 
a CP-odd particle, its mass has to be larger than at least 3m n for the decay into final states 
without electromagnetic particles. Therefore, if m a < Sm^, a — > 77 is expected to be the 
dominant decay mode. In this case, we can use the Eq. (|3.20|) again with q;sm = ^em- In 



the case 3m, < m a ^ 1 GeV, estimation of the decay rate is quite difficult, and we will not 
discuss this case further in this paper. 

In Figs. f| and ||, we plotted fl^h 2 on the v vs. Bq plane with O = and N 5 = 2. In 
the calculation, we used Eq. ( |4.17|) with ^sin^o = 1, and assumed that a — >• 77 and a —> gg 
are the dominant decay modes of the i?-axion for m a < "im^ and m a > 1 GeV, respectively. 
In these figures, we shaded the region where the reheating temperature becomes lower than 
1 MeV. We also show the contours of the constant v/F^/ 2 , which has to be larger than 
~ 1 — 10. 

Let us first discuss the case of m a < 3m n , where the i?-axion decays into the photon 
pair. This is the case for the smaller value of v (i.e., 10 8 — 10 9 GeV). If the i?-axion 
decays into the photon pair, the reheating temperature becomes relatively low. However, 
if Bq is large enough, there is still a parameter region where the reheating temperature is 
high enough (T R > 1 MeV). At the same time, can be of O(0. 1) or smaller, and hence 
the moduli field can be diluted enough by the decay of the i?-axion. 
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Figure 4: Constant Vt^ti 2 contours on the v vs. Bq plane with 0o — M*, N$ = 2 and 
k Q = 1. Region with Tr < 1 MeV is lightly shaded, and darkly shaded region corresponds 
to 3/77.^ < m a < 1 GeV. Contours of constant v/F^ are also shown in dotted lines (1, 10, 
and 100, from left to right). 
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Figure 5: Same as Fig. B, except for k = 3. 
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One problem of this case may be the non-perturbative effect discussed in Section ^. In 
the case of v < 10 9 GeV, non-perturbative effect becomes sizable in the Kahler potential, 
and hence the potential of X cannot be well understood. Therefore, we need to make 
a dynamical assumption so that the SUSY breaking vacuum exists. However, since the 
non-perturbative effect respects the i?-symmetry, important properties of the .R-axion are 
unchanged. Therefore, once we assume the existence of the SUSY breaking minimum, our 
arguments are unchanged. One may also worry about the stability of the vacuum if there 
is a SUSY preserving true minimum at the origin. However, T R > 1 MeV can be realized 

1/2 

even with v/F^ £ 10. Furthermore, this problem can be evaded in models without SUSY 
preserving minimum. 

For a larger value of v (i.e., v ^ 10 9 — fO 10 GeV), the i?-axion becomes heavier than 
~ 1 GeV, and decays into the gluon pair. As we can see in Figs. [| and ||, if the i?-axion 
decays into the gluons, the density parameter of the moduli field becomes relatively larger. 
This is because the decay rate is more enhanced, resulting in higher reheating temperature. 
Even in this case, however, ~ 10 is possible for O ~ M*. Reduction of by a factor of 
about 10 may not be a serious problem. For example, suppression of the initial amplitude 
of the moduli field by a factor of 3 or so is enough for < 1, since is proportional 
to 0o- We can naturally imagine an accidental cancellation of this level. Notice that the 
non-perturbative effect on the Kahler potential is not important in (most of) this case. 



4.4 Effect of the Real Part 

Next, we discuss the effect of the decay of the real part a. For this purpose, it is convenient 
to use the relation m e f[(p 2 oc R~ 3 (with ip = X and 4>). Comparing this quantity at \X\ ~ Xo 
and \X\ ~ v , we obtain 



m 



3/20 2 




(4.20) 



The important point is that the effective mass for the moduli field is always constant of 
0(7713/2), while that for X varies from ~ m 3 / 2 to m a which is of the order of the SSM scale. 

Once X gets trapped in the SUSY breaking minimum, a obeys the parabolic potential, 
and hence Eq. ( f4.20| ) leads 



Pa 




(4.21) 



where p a is the energy density of a. 

When the expansion rate becomes comparable to T a , a decays. As discussed in Section |[ 
a dominantly decays into the .R-axion pair. The number density of a, n a , is related to the 
energy density as p a = m a n a . Therefore, by using Eq. Q4.21]) with \Xq\ ~ M#, we obtain 



^dcc 
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where n^ ec is the number density of the i?-axion produced by the decay of a. Notice that 
this ratio is constant even with entropy production. 

As the Universe expands, the emitted .R-axions are red-shifted and eventually become 
non-relativistic. Then, when they decay, they also contribute to the entropy production to 
dilute the moduli field. Comparing Eq. ( |4.22|) with Eq. ( |4.6|) , we see that the number of the 
i?-axion produced by the decay of a is comparable to that of the coherent mode. Therefore, 
the dilution by this incoherent i?-axion is of the same order of that by the coherent mode, 
and the results given in the previous subsection are almost unchanged. 

Finally, we discuss the effect of the decay mode into the gravitino pair. As discussed in 
Section [| a and a may decay into the gravitino. Since the gravitino is stable in GMSB, the 
mass density of the gravitino, P3/2, should not exceed the closure limit. 

Let us consider the gravitino production due to the decay of a as an example. In order 
to discuss the mass density of the gravitino by the decay of a (and of the i?-axion), it is 
convenient to consider the ratio m 3 / 2 n 3 / 2 / p^ ~ P3/2/P0, where n 3 / 2 is the number density 
of the gravitino. This ratio can be easily calculated from Eq. ( |4.21| ). When o decays, the 
number density of the gravitino is given by 

n 3/2 ~ Br(cr WVv) x n * ~ Br ( a VV^V) x — • ( 4 - 23 ) 
Combining this equation with Eq. ( [4.21 ) and \Xq\ ~ M*, we obtain 



Importantly, this ratio remains constant even after a large entropy production. Once the 
gravitinos are red-shifted to be non-relativistic, the above ratio becomes P3/2/ 1 p<j>- As a result, 
the density parameter of the gravitino, fi 3 / 2 , is related to as 

n 3/2 ~ Br(cr - ^) x % (Jpj . (4.25) 

It is also straightforward to check that the density parameter of the gravitino from the decay 
of the i?-axion is given by a similar formula (with extra coefficient of O (1)) with the relevant 
branching ratio for the i?-axion. 

With the branching ratio given in Eqs. ( |3.25| ) and ( |3.26| ), we can see that the mass density 



of the gravitino is small enough. By using the formula for m^, Br(<7 — > ijj^ijj^ is estimated 
to be of 0(1CT 4 ). Furthermore, in most of the parameter region we are interested in, the 
i?-axion mass is smaller than about 10 GeV, and hence Br(a — > VvVv) ^ s a ^ m ost of O(10 -2 ). 
Therefore, even with fi^ x (0 O /M»,)- 2 ~ 10 for v > 10 9 - 10 10 GeV (i.e., for m a > 1 GeV), 
the energy density of the gravitino is small enough. 
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4.5 Baryogenesis and the Flat Direction in SSM 

So far, we have discussed the dilution of the moduli field by the decay of the .R-axion. In 
this scenario, one may worry about the baryogenesis, since the reheating temperature seems 
to be too low to generate the baryon asymmetry. In this subsection, we briefly comment 
that enough baryon asymmetry can be generated by the Affleck-Dine mechanism (For 
a detailed discussion, see Ref. 0.) 

Important assumptions for the Affleck-Dine baryogenesis are that SSM flat direction 
(Pssm (for Affleck-Dine baryogenesis, we call it Affleck-Dine field) has a large amplitude in 
the early Universe, and that there is baryon-number breaking term in the potential of y?ssM- 
With these assumptions, non- vanishing baryon-number can be generated as the Affleck-Dine 
field evolves. Physics in this mechanism is basically the same as that in the case of the 
.R-number generation. ^ 11 ! Since the potential for the Affleck-Dine field is dominated by the 
supergravity contribution for y?ssM ~ t ne ratio of the energy density of the moduli field 
to the baryon number density jib is estimated as [f| 



n B 




(4.26) 



where we applied a similar argument as in the i?-axion case. Here, we assumed that the 
initial amplitude of y^ssM is of 0(M*) and that source of CP violation is of 0(1), which are 
corresponding to £ ~ 1 and sm6 ~ 1 in the i?-axion case, respectively. Notice that this 
ratio is constant after the moduli and Affleck-Dine fields start to move. 

With Eq. ( |4.26j ), we can estimate the baryon number density of the present Universe, 
once we fix the current mass density of the moduli field. As we discussed in the previous 
subsections, present mass density of the moduli field depends on the magnitude of the entropy 
production. In this subsection, we just assume the moduli field is diluted enough by the decay 
of the i?-axion. Then, adopting ^ 1, baryon-to-entropy ratio is estimated as [0] 

^ < 4 x 1CT 5 x h-> ( m f) " ( # V * . (4.27) 



Comparing the above relation with the baryon number density required from the big-bang 
nucleosynthesis {ti-q/s ~ 0(lCU n ) f29|), we can see that enough baryon asymmetry can 



remain even if there is a large entropy production to dilute the moduli field. In fact, as 
can be seen in Eq. ( [4.27|) , baryon-to-entropy ratio may be too large for some value of the 
gravitino mass, if we adopt a naive initial condition for the Affleck-Dine field. However, 
this problem may be solved by adopting a smaller value of the source of CP violation, or a 
smaller value of the initial amplitude of the Affleck-Dine field. Therefore, we do not worry 
about this issue. 

Even apart from the baryogenesis, effect of the SSM flat direction may be interesting, 
since it may produce a large entropy. The potential for the flat direction has a similarity to 

# n The structures of the baryon-number breaking terms may be different. However, it does not affect the 
following discussions, as far as the initial amplitude of </?ssm is of 0(Af„). 
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that of the a field; it is parabolic around the minimum (i.e., around the origin), then the 
potential is lifted only logarithmically once the amplitude becomes larger than the messenger 
scale, and finally, supergravity effect dominates when the amplitude is close M*. Therefore, 
its effective mass varies from ~ ^3/2 to the SSM scale, if <^ssm changes its amplitude from 
~ M* to ~ 0. By applying the similar argument as in the a case, we can estimate the density 
parameter of the moduli field: 

n *" ~ 2 x 10 x IPSvj (itevJ [Wkv) [w,) • (428) 

where T SSM is the reheating temperature due to the decay of the flat direction. Comparing 
this result with Eq. ( |2.1| ), we can see that a large entropy production is possible even from 
the flat direction in SSM. 

However, usually, the reheating temperature is of O(10 TeV), and hence the density 
parameter becomes typically of O(10 3 ) even for m 3 / 2 ~ 1 keV. Therefore, it seems difficult 
to dilute the moduli density enough in a simple scenario, unless we come up with a model 
with very light gravitino. However, if the reheating temperature can be somehow lowered 
down to ~ 10 GeV, dilution due to the decay of the SSM flat direction may be enough to 
solve the cosmological moduli problem. For example, Pauli blocking may delay the decay 
of the flat direction, as pointed out in Ref. fl3~0"| , although the baryon asymmetry may be 
generated too much if we naively apply the argument in Ref. pO . 



5 Discussion 

In this paper, we have discussed the cosmology based on the direct gauge mediation model 
with the inverted hierarchy mechanism. In particular, we have studied the implication of 
the SUSY breaking field on the cosmology of GMSB. 

If the SUSY breaking flat direction X initially has a very large amplitude of 0(M»), it 
can be a source of the large entropy production. In particular, once the amplitude of X 
becomes a few orders of magnitude smaller than M#, potential for X is dominated by the 
logarithmic piece. Then, the energy density of X decreases more slowly than those of the 
scalars with quadratic potential, and the SUSY breaking field may play a very important 
role in cosmology. 

In particular, the entropy production by the decay of the .R-axion may be so large that 
the energy density of the moduli field can be diluted enough. Therefore, the direct gauge 
mediation models with the inverted hierarchy mechanism contain a natural candidate of the 
large entropy production to solve the serious cosmological moduli problem. We have also 
seen that enough baryon asymmetry can be generated by the Affleck-Dine mechanism even 
with this large entropy production. 

In our discussion, we mainly focused on the case of the direct gauge mediation model 
with the inverted hierarchy mechanism. However, the scenario discussed in this paper can 
be applied to a larger class of models, since the most important building block is just the 
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logarithmically lifted potential of the SUSY breaking field for \X\ ^ t>. Consequently, if a 
direct gauge mediation model uses the mechanism of SUSY breaking proposed by Izawa 
and Yanagida, and by Intriligator and Thomas |14[], it automatically contains a reasonable 
source of large entropy production. Of course, it is non-trivial for such a model to stabilize 
the potential of the SUSY breaking field. The inverted hierarchy mechanism provides one 
attractive mechanism for the stabilization [0, §, [|, |10j- Other approach may be to use a 
non-perturbative effect on the Kahler potential (Tl[ . 



One may worry about the astrophysical constraints on the .R-axion, as in the case of the 
QCD axion. However, since the .R-axion has a larger mass than the QCD axion, constraints 
on the .R-axion are much weaker. Constraints from the cooling of the horizontal-branch (HB) 
stars [31] are evaded, since we consider .R-axion heavier than 100 keV (see Eq. (|3.19| )) while 



the core temperature the HB stars are typically of O(10 keV). Furthermore, the .R-axion 
does not affect the background UV light p2| , since it has already decayed away. Constraints 
from SN1987A are more non-trivial. In order not to affect the cooling of SN1987A, QCD 
axion with the decay constant from ~ 10 6 GeV to ~ 10 9 GeV is forbidden ||33|| . However, 
for this range of v , the .R-axion mass is (almost) always larger than the core temperature of 
SN1987A (i.e., O(10 MeV)), and hence the emission of the .R-axion is suppressed. (However, 
small parameter space around v ~ 10 6 GeV and Bq ~ 10 4 GeV may be excluded, though 
the reheating temperature is too low in this region. ) When v < 10 6 GeV, the R- axion is 
thermalized enough in SN1987A, and it does not affect the cooling process. The QCD axion 
with small decay constant may be detected in water Cerenkov detectors EM. However, the 



.R-axion cannot be constrained with this method, since .R-axion decays before reaching the 
earth. This fact suggests another constraint on the .R-axion; if the emitted .R-axions decay 
into the photons on the way to the earth, apparent luminosity of SN1987A may be increased. 
Therefore, for our scenario, light .R-axion is potentially dangerous. However, the estimation 
of the .R-axion flux is very complicated, in particular since the .R-axion may have a mass 
comparable to the core temperature of SN1987A. Therefore, it is an open question which 
parameter region is excluded from this argument. Notice that, if the R- axion mass is heavier 
than of O(10 MeV), this problem can be evaded thanks to the Boltzmann suppression. 

In this paper, we have not paid attention to the primordial inflation, since it is beyond 
the scope of this paper. Of course, it is important to find a viable candidate of the inflaton 
for the primordial inflation, and some efforts are made in this issue [35|. Here, we just 
mention that, in our scenario, inflaton for the primordial inflation does not have to decay 
into the particles in the SSM sector. Since the background radiation and baryons in the 
present Universe originate to the decay of the SUSY breaking field X (and probably, to the 
decay of the Affleck-Dine field), primordial inflation is not required to reheat the SSM sector. 
In an extreme case, inflaton may decay only into particles in the hidden sector. Even if the 
energy density of the Universe is once dominated by that of the hidden sector particles, it is 
eventually diluted by the entropy production by the decay of the SUSY breaking field. This 
fact relaxes the conventional requirements on the inflaton which is usually required to decay 
into the SSM particles. 

It is interesting to consider candidates of the CDM in this scenario. Because of the low 
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reheating temperature after the decay of the .R-axion, thermal productions of any known 
candidates are inefficient. In this case, a possible candidate is the coherent oscillation of the 
moduli field. Indeed, if the energy density of the moduli field can be diluted to be fi^ ~ 1 
by the decay of the i?-axion, they can be a viable candidate of the CDM. This scenario is 
constrained by the line spectrum of the background cosmic X-ray emitted from the decay 
of the moduli field [36[. Even if the lifetime is much longer than the age of the Universe, 



some fraction of the moduli field has already decayed into the photons, and it contributes to 
the background X-ray spectrum. Due to the negative observation of such a line spectrum, 
moduli field heavier than about 200 keV is forbidden, if Q^h 2 = 1 p6| . | ff x j In the direct 
gauge mediation model, the gravitino mass (i.e., the moduli mass) can be lighter than about 
100 keV, and hence the moduli can be a good candidate of the CDM. 

Cosmological implication of the i?-axion in other classes of GMSB is another interesting 
issue. In general, dynamical SUSY breaking requires spontaneously broken i?-symmetry [[37] . 
Therefore, all the dynamical SUSY breaking models contain i?-axion in the low energy 
spectrum. Then, there is a possibility of large i?-number generation by the mechanism we 
discussed, if the SUSY breaking field has a large initial amplitude. However, the interaction 
of the i?-axion is model-dependent, and in some case, its decay rate may be much more 
suppressed. In this case, it causes a cosmological difficulty, since the reheating temperature 
after the decay of the i?-axion becomes too low for the big-bang nucleosynthesis. Of course, 
this problem itself can be always evaded by assuming a small initial amplitude of the SUSY 
breaking field. 

Since different models introduce different sets of new particles which have various prop- 
erties, detailed cosmological scenario is model-dependent. Therefore, one should always keep 
in mind that the SUSY breaking field (and all the new degrees of freedom) may change the 
conventional arguments on the cosmology based on supersymmetric models. In some case, it 
may cause a serious cosmological disaster, but in other case, as we have seen, it may provide 
a natural and well-motivated solution to several serious cosmological difficulties. 
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# 12 However, if the interaction of the moduli field with the photon has extra suppression, lifetime of the 
moduli field gets longer and this constraint becomes less stringent. In this case, heavier moduli field is still 
allowed. For example, this may happen if the moduli-photon-photon vertex is induced by loop effects. 
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A Example of the Model 



In this Appendix, we show an example of the direct gauge mediation model in which vacuum 
is stabilized by the inverted hierarchy mechanism. The model, which is originally proposed 
in Ref. ||, is based on the symmetry SU(2) B i x SU(2) B 2 x SU(2) S x G S m- In this model, 
SU(2)bi is a gauge group which stabilize the minimum of the potential, while SU(2) S provides 
a strong gauge interaction which induces gaugino condensation to break SUSY. Furthermore, 
Gsm is the standard model gauge group. The particle content is shown in the Table |I[ 
Superpotential in this model is given by 

W = Vq^QQ + y 3 T,q s q 3 + y 2 T,q 2 q2 + Vi^qiQi- (A.1) 

With this superpotential, we concentrate on the flat direction of X ~ (detE) 1 / 2 . Along X, 
we parametrize 

S = i=diag(X,X). (A.2) 

Once E gets this VEV, Q and Q acquire a mass of tuq ~ A=ygX. For jj, ^ iuq, SU(2)s is a 
pure SUSY Yang-Milles theory, and gaugino condensation induces the superpotential of the 
form W e fi = 2Ag ff , where A eff is the strong scale of SU(2) S below the mass scale of Q and Q. 
By matching the strong scales for the theory below and above tuq, we obtain 

W eS = V2y Q XA 2 , (A.3) 

where A is the strong scale for the theory above tuq. Since the superpotential is linear in 
X, F-component of X has a VEV of F x = ^/2ynA 2 , and SUSY is broken. 



The minimum of the potential is determined by the inverted hierarchy mechanism [[13 
At the tree level, potential for X is completely flat, and hence the scale dependence of the 
wave function normalization of E determines the position of the minimum. In this case, the 
potential for X is given by 

F 2 

V = — r, (A.4) 

where Z^,(X*, X) is the wave function normalization of E which is evaluated at \i — \X\. At 
the 1-loop level, RGE for is given by 

rflnZs 1 /3 2 3 



dt 16tt 



2^bi + - %J - 3 % 2 - 2y 2 2 - y 2 ) , (A.5) 



where g B i and gB2 are the gauge coupling constants for SU(2) B i and SU(2) B2 , respectively. 
Thus, X has an extremum at X — v, where 

+ 9I2) = + 32/3 + 2y 2 2 + yl (A.6) 
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Table 1: Particle content of a direct gauge mediation model given in Ref. |§. 



In order to see whether this is a minimum or a maximum, it is convenient to estimate the 



mass of the real part of X, which we call q.| yl3 | For simplicity, we consider the case where 



9B2, 1/3, and y<i are small. (Even in the general case, the following discussion is qualitatively 
correct.) Then, mass of a is given by 



.2 ^ ' 4 | n/| 4 0J( _2 „.2 e„2„.2\ 



, 1C)7r2 )2 v 4 + 24yJ - 24^ - 6^ j -A, (A.7) 



where g$ is the gauge coupling constant of SU(2)s. Importantly, gs is usually large in order 
to induces the strong dynamics to break SUSY. Then, for large yq, g$yq term becomes so 
large that m 2 a becomes negative. (For large enough yq, m\ may become positive, but yq 
blows up below the Planck scale.) Numerically, yq cannot be larger than 0.2 - 0.3 for the 
positivity of m 2 a . For v ^ 10 9 GeV, solution to Eq. (|A.6|) with positive m\ can be found with 
reasonable values of the coupling constants. 

For v < 10 9 GeV, we cannot neglect the non-perturbative effects, as we discussed in 
Section || In this case, Kahler potential is dominated by the non-perturbative piece, and 
it is unclear whether there can be a minimum. However, since the non-perturbative effects 
are not well understood, there is a possibility to have a stable minimum even with the non- 
perturbative effect. Therefore, for v ^5 10 9 GeV, we make a dynamical assumption so that 
the stable SUSY breaking vacuum exists. Notice that, in this case, upper bound on yq is 
irrelevant. 

Once the SUSY is broken and the VEV of X is fixed, SUSY breaking is mediated down 
to the SSM sector by integrating out messengers, and g« {i = 3,2). Since these are the 
only superfields with standard model quantum numbers which couple to the SUSY breaking 
field, SUSY breaking masses obey the well-known mass formula ||38|| . Notice that, in this 
model, N 5 = 2 (see Eq. (KW)). 



# 13 Imaginary part is pseudo-Nambu-Goldstone boson, and its mass is from the supergravity effect, as 
discussed in Section [| 
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Several remarks are in order. First, for fM^v, diagonal SU(2) remains if SU(2)b2 is 
gauged. For our cosmological scenario, this may be dangerous, since the i?-axion also couples 
to this gauge field. If .R-axion dominantly decays into this gauge field, energy density of 
the standard model particles cannot be generated enough since the gauge field for this 
diagonal SU(2) does not couple to the SSM sector. This problem can be evaded if SU(2) B 2 
interaction is weak enough (i.e., much weaker than the electromagnetic interaction, or maybe 
not gauged). Another comment is on the true vacuum of the potential. In fact, this model 
has a SUSY preserving true vacuum at the origin (\X\ = 0). However, tunneling rate from 
the SUSY breaking minimum to the true minimum is small enough, if the SUSY breaking 
vacuum is far enough away from the confinement scale A; numerically, v/A^10 ||. 



B Shift of the Minimum of the Moduli Potential 

In this Appendix, we consider how the minimum of the moduli potential shifts when the 
SUSY breaking field X is oscillating. When the Universe is dominated by X, extra terms are 
induced in the moduli potential by supergravity effects, and they may change the minimum 
of the moduli potential. If this shift is too large, it may change our argument, since the 
moduli field may oscillate around a shifted minimum. In this Appendix, we see that this 
effect is not significant for our case, and that our naive calculations are relevant. 

First, let us consider possible modifications of the moduli potential in the presence of 
X. In supergravity, there can be two effects. One is from the non- vanishing VEV of X; 
since X and may have Planck-suppressed interactions, potential may have terms which 
are proportional to the powers of (\X\/M*). The other is from the expansion rate H induced 
by the condensation of X; since the scalar potential contains a term of the form ~ e K / M * V, 
non- vanishing potential energy induces terms proportional to H 2 [p"5|j . 

With these effects, linear term is induced in the moduli potential, which shifts the mini- 
mum of the potential: 

V{(j>) ~ ml /2 (j>*4> - ml/S^ + h-c), (B.l) 

where we define the origin of the moduli field so that VEV of <p vanishes for the empty 
background. Here, the second term in Eq. (|B.1|) is the induced term, and 4> is given by 



~ max 



\X\, {H 2 /mi /2 )M* 



fB.2) 



(In this section, we neglect 0(1) coefficients which do not change our argument.) Notice 
that 4> is the shifted minimum of the potential, and that it is time-dependent With the 

# 14 Other terms (higher order terms) are less significant for our argument, and they do not change the 
following discussion. For example, if the potential has a term of the form H 2 (f> 2 , solution to the equation 
of motion contains a term of 0(H 2 /m 2 / 2 )(j) sc, with <j> osc being the solution to the equation of motion with 
H = 0. However, this is much smaller than the original amplitude 4> OS c since H <C TI3/2, and hence this 
effect is negligible. 
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above potential, equation of motion for is given by 



+ 3F0 + m^ /2 (0-0) = 0. (B.3) 
Solution to this equation can be written as 

+ 50, (B.4) 



y osc 



where OSC is the oscillating solution with = 0, while £0 is a perturbation induced by the 
new terms. Notice that OSC obeys the original equation of motion: 

OSC + 3H<j) osc + ml /2 (j) osc = 0, (B.5) 

and hence its averaged amplitude is proportional to R~ 3 ^ 2 . Thus, OSC obeys the behavior 
discussed in Section |. 

If ~ |X|, the shift cannot be larger than the averaged amplitude of OSC (see Eq. ( |4.11]) ). 
In this case, the original amplitude is always larger than the shift of the minimum, and hence 
the extra contribution is negligible. Thus, in the following discussion, we concentrate on the 
case where is dominated by the Hubble-induced term. 

In order to consider the Hubble-induced term, we approximate the potential of X as 

v ^(rlr m ' /2M '( ln ^?) 2 ' (R6) 

with (2 being a constant. This potential has a minimum at |X| = v, and increases loga- 
rithmically for large |X|. Therefore, this potential reproduces the important feature of the 
potential of X (especially for |X| ~ v). When the Universe is dominated by X, expansion 
rate is estimated as 

H ~ 7^ Lm 3/2> ( B - 7 ) 



with 



ln^. (B.8) 



v 2 

Importantly, this expansion rate is much smaller than the gravitino mass. With the above 
expansion rate, we denote 

0=^M„ (B.9) 

m 3/2 



where kn is an unknown constant expected to be of 0(1). We consider the case where the 
energy density of the Universe is dominated by that of X, so H decreases as the amplitude 
of X approaches to v. 
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As shown in Eq. 



, solution to Eq. 



is given by the sum of the oscillating solution 



OSC which is from the non-perturbed equation of motion and 50 induced by extra terms with 
non-vanishing H. We have already understood the behavior of OSC , so now we consider 5(f>. 

The important point in deriving 8(f) is that the expansion rate H is much smaller than 
the gravitino mass m 3 / 2 - Because of this, we can expand <50 by powers of (H 2 /m^ 2 ). First, 

let us consider a simple case where H is proportional to H 2 : 



H 



-c H H 2 



with cjj being a constant of 0(1). For example, if X has a parabolic potential, c# 
In this case, solution to Eq. ( [B.3| ) is obtained as 



+ 



k H H 2 



rn 



- Qc H (c H - r 



k H H 4 



+ 



H 6 



3/2 



m 



3/2 



m 



3/2 



(B.10) 
= 3/2. 



fB.ll^ 



Notice that the term of 0(H 2 /m^ 2 )M Jf is exactly equal to (f>. In the case where the potential 
of X is logarithmic like Eq. ( P-6|) , formula for H is slightly different: 

1 



H 



-H 



In this case, is given as 



+ 



k H H 2 36(5 + 2L - L 2 ) k H H 4 



in 



3/2 



o 



rn 



3/2 



m 



3/2 



(B.12) 



(B.13) 



and the 0(H 2 /m^ 2 )M* term agrees with again. In general, as far as H m 3 / 2 and 



H ^ 0(H 2 ), the leading correction is always equal to (f>. This is because, in Eq. flB.3|) , first 



two terms become of 0(H 4 /m^^M* 



and hence (f> has to be cancelled out by 0(H 2 /m^i 2 )M t 



4 )M*. 

smoothly goes to as the Universe expands, (f> term in (f> is harmless. In other 



As a result, the deviation from the shifted minimum is always of 0(H /rn 3 / 2 



term in q 
Since 

words, in the early stage, <f) oscillates around the shifted minimum 0, but this minimum 
approaches to the true minimum as H — > 0. Therefore, we just have to consider the deviation 
from = 0. 

In our situation, potential of X changes its behavior at |X| ~ v ; it is logarithmic for 
\X\ ZjV, and parabolic potential is relevant once X is trapped in the SUSY breaking mini- 
mum. Thus, the solution to Eq. (p.3j ) changes its behavior at \X\ ~ v. For example, if we 
match two cases at L — 1 (though the matching point is quite uncertain), OSC is shifted as 



45 k H H 4 



-M* 



m 



(B.14) 



3/2 



Notice that this shift is not 0(H 2 /m|/ 2 )M*, but 0(H 4 /m\, 2 )M*. In general, shift of this 
order may be possible, especially when the potential changes its behavior. However, shift 
cannot be 0(H 2 /m\, 2 )M^. 
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If this shift is larger than the averaged amplitude of O sc, our argument may break down. 
Therefore, we require 



W>4-M*. (B.15) 
m %2 



When \X\ ^v, OSC is proportional to IX) 1 / 2 while H depends on \X\ only logarithmically. 
On the other hand, once X is trapped in the SUSY breaking minimum, OSC and H are both 
proportional to (|X| —v). Therefore, this constraint becomes most stringent when |X| ~ v. 
Since the expansion rate is much smaller than the gravitino mass (see Eq. (|B.7|) ), above 



constraint is very weak. Numerically, this requires 100 GeV, even for £ 2 ~ 1. Of course, 
if we adopt smaller value of (2, lower bound on v becomes less severe. Therefore, in our case, 
this constraint is weak enough, and Hubble-induced term does not change our argument. 



C Evolution of tir 

In this Appendix, we discuss how the i?-number density evolves with time. In particular, 
since i?-symmetry breaking terms exist in the potential, we need to know the evolution of 
tir with them. 

We consider a case where X oscillates with a large amplitude. In this case, frequency 
of the oscillation is roughly given by ~ m e s, where the effective mass m e g is given by (see 
Appendix 0) 

2 _ 1 ( 1 9V 1 dV \ 
meff "2 \X*dX + XdX~*)- 

Notice that m e fj ^> H when X is oscillating. Therefore, we consider the time scale <^ 
5t H 1 , for which we neglect the change of \X\ and H. With the expansion of the 
Universe, the amplitude of X decreases as 

\X\R P = const., (C.2) 

with p being a positive 0(1) constant. (For example, p = 3/2 for parabolic potential, and 
p = 3 for logarithmic one.) Therefore, we approximate the motion of X as 

X ~ \X \e {irncS - pH)t . (C.3) 

With the i?-symmetry breaking potential 

V n ~ + x) x f( x * x / M * )' ( c - 4 ) 

with f(x) = k + k\x H — equation for the evolution of ur = i(X*X — X*X) is given by 

h R + 3Hn R = 2Fllm(X/M,) x /(|X| 2 /M 2 ). (C.5) 
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For the following discussion, it is more convenient to consider the i?-number in a comoving 
volume. For this quantity, the above equation leads to 

= R 3 x 2F*1m(X/M m ) x f(\Xf/M*). (C.6) 

We need to solve the above equation to obtain the resultant i?-number. For this purpose, 
we first take the average of the right-hand side of this equation for the time scale <C 
St <C H~ x . For this time scale, we approximate \X\ and H to be (almost) constant. On the 
other hand, the change of ImX is extremely non-adiabatic. In the flat background (i.e., if 
H = 0), average of ImX is supposed to vanish since X is in a periodic motion. However, in 
the actual situation, H is non- vanishing. By using Eq. flC3[ ), average of ImX is estimated 
as 

(ImX)~\X\— (C.7) 

m eS m eS dt 

where we related the expansion rate H to d\X\/dt by using Eq. (|C.2|) . In Eq. ( |C.7| ) and here- 
after, we neglect possible 0(1) coefficients since they do not change the following argument. 
Combining this equation with Eq. ( |C.6| ), we obtain 

„ fl 3 x ^Jl iWMl . (c.8) 

d\X\ M* m e s 

By using the fact that m e s\X\ 2 R 3 is a constant of motion, we integrate the above equation 
from X = Xi to X = X f (with X { > X f ): 

R 3 n R ~ m eS \X\ 2 R 3 x x y J^-( |A / - - '7' - . (C.9) 

If the second term in the parenthesis in Eq. flU.9 ) wins the first term, integration in large 
\X\ region is more important, and major part of the i?-number is generated when X starts 
to move. On the other hand, if the first term is dominant, we cannot neglect the i?-number 
generation in the later stage. 

If the potential of X is dominated by the logarithmic piece, m e g is proportional to 
In this case, i?-number generation at large amplitude is more important for n > 0. As a 
result, even if there is a linear i?-symmetry breaking term in the potential, R- number is 
generated when X starts to move, and i?-symmetry is conserved with a good accuracy for 
a smaller value of |X|. On the other hand, for parabolic potential, m e a is a constant. If 
the R- number violating potential is dominated by the linear term (n — 0), contribution at 
small \X\ becomes important. However, in our scenario, we assume that the linear term is 
suppressed enough when X starts to move, and that the i?-number violating effect starts 
with cubic term (n = 1). In this case, i?-number asymmetry is again generated when X 
starts to move. 

In the actual situation, X starts to move with a quadratic potential, and at some stage, 
logarithmic piece takes over. With the assumption that the linear term is suppressed enough, 
i?-number is generated when X starts to move, and afterwards, i?-number in the comoving 
volume is conserved well. 
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D Scalar Field in the Expanding Universe 

In this Appendix, we derive a convenient formula for the evolution of the scalar field ip in 
periodic motion. For simplicity, we consider the case where the amplitude of the scalar 
field is (almost) constant in a time scale of the periodic motion and also the potential for ip 
depends only on \<p\. 

From the virial theorem, we obtain 

where /C = ip*ip is the kinetic energy of ip, and the bracket represents the time average. 
The field equation for ip is given by 

dV 

^ + 3^+^ = 0. (D.2) 

Multiplying this equation by ip* , and using the definition of /C, we obtain 

/ dV dV\ 

Now, we are at the position to consider the evolution of the scalar field. For this purpose, 
we define 

(D - 4) 

and consider the evolution of (S 2 ). By taking the derivative of (S 2 ) with respect to time, 
we obtain 

*W = ldW?( OV 0V_\ L f 3V 0V\\ 

dt \ dt \jd<p * d v * ) m dt yd? * dip* J I 

ldW? ( dV *dV\ 2 -\ 
= -12/j(|v?| 2 /C 

= -6H(S 2 ), (D.5) 

where we used the fact that the potential V is a function of \<p\. By solving the above 
equation, we obtain 

(S 2 }R 6 = const. (D.6) 
The scalar field evolves by following this relation. 
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For a more intuitive understanding, it is convenient to define the "effective mass" from 
the potential V: 

2 f / f dV 1 dV \ 

With this effective mass, evolution of the scalar field is given by 

m e ff|<^| 2 .R 3 = const. (D.8) 

For example, in the case of parabolic potential V = m^\ip\ 2 , m eS does not depend on ip, and 
hence \<p\ scales as i?~ 3//2 , while \(p\ oc R^ 1 for quartic potential V oc |</9| 4 . 
Notice that, in the flat space (H = 0), 

y = cp e ±im ^\ (D.9) 

satisfies the equation of motion of ip for any value of ipo, if the potential of ip depends only 
on \<p\. (In Eq. (|L).9|) , m e fj is evaluated at ip = <po.) Therefore, m e s can be understood as a 
frequency of the periodic motion. 

Since m e g|<y9| 2 is proportional to the volume of the phase space for the periodic motion, 
a physical interpretation of Eq. ( |D.8| ) is that the phase space volume in a comoving volume 
is conserved as the Universe expands. 
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